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Abstract 



o 

^5 ' In this article, we assume that there exists a scalar hidden charm tetraquark state 

in the t^'^XcI invariant mass distribution, and study its mass using the QCD sum 



S ' rules. The numerical result Mz = (4.36 ± 0.18) GeV is consistent with the mass of 

the ^(4250). The Z(4250) may be a tetraquark state, other possibilities, such as a 



hadro-charmonium resonance and a D^D^ + D^ D\ molecular state are not excluded. 
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^' 1 Introduction 

Recently the Belle collaboration reported the first observation of two resonance-like struc- 

\^ , tures in the 7r"'"Xci invariant mass distribution near 4.1 GeV in the exclusive B^ -^ K~7r~^Xci 

^ I decays [1]. The Breit-Wigner masses and the widths are about Mi = 4051 ± lA'^fl MeV, 

Ti = 821^^1^^ MeV, Ma = 42481^^ +^f MeV and Fa = 177+|^l^f MeV (thereafter we 

lO ■ will denote them as Z(4050) and ^(4250) respectively). The significance of each of the 

''^ . TT'^'Xci structures exceeds 5a, including the effects of systematics from various fit models. 

r^ I Their quark contents must be some special combinations of ccud, just like the Z(4430), 

QQ ■ they cannot be the conventional mesons [2j. 

(^ . The Z (denote the Z(4050) and Z(4250)) may be loosely deuteron-like bound states 

(molecules) of the charm mesons or compact nucleon-like bound states of the diquark- 
antiquark pair. The spins of the Z are not determined yet, they can be scalar or vector 
states. 

5^ \ In the meson-exchange model, the Z(4050) is probably a loosely molecular state D*D* 

with J^ = 0^ [3j and the Z(4250) is unlikely an S'-wave DiD or DqD* molecular state 
[4|, while the SU{3) chiral quark model indicates that the Z(4050) is unlikely an 5- 
wave D*D* molecular state [5]. In Ref.[6l[7j, the authors study the mesons Z(4050) and 
Z(4250) as the D*D* and D+D^ + D+D^ molecular states with J^ = 0+ and J^ = 1" 
respectively using the QCD sum rules, and draw the conclusion that the D*D* state 
is probably a virtual state which is not related with the Z(4050) and the Z(4250) is a 
possible DfD^ + D^D^ molecular state. In a relativistic quark model, the Z(4250) can 
be tentatively interpreted as the charged P-wave 1~ tetraquark state SS or as the P-wave 
0~ tetraquark state {SA it .Syl)/v2 [8], where the S and A denote the scalar and axial 
vector diquarks respectively. 

The colored objects (diquarks) in a confining potential can result in a copious spectrum, 
there maybe exist a series of orbital angular momentum excitations; while the colorless 
objects (mesons) bound by a short range potential (through meson-exchange) should have 
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a very limited spectrum. In the heavy quark hmit, the c quark can be viewed as a static 
well potential, and binds the light quark q to form a diquark in the color antitriplet channel. 
We take the diquarks as the basic constituents following Jaffe and Wilczek [9l [10]. The 
heavy tetraquark system could be described by a double-well potential with two light 
quarks q'q lying in the two wells respectively. 

In Refs.pn [T2| [T3l [Til [15], Maiani et al take the diquarks as the basic constituents, 
examine the rich spectrum of the diquark-antidiquark states from the constituent diquark 
masses and the spin-spin interactions, and try to accommodate some of the newly observed 
charmonium-like resonances not fitting a pure cc assignment. The predictions depend 
heavily on the assumption that the light scalar mesons ao(980) and /o(980) are tetraquark 
states, the basic parameters (constituent diquark masses) are estimated thereafter. 

In Ref . [I6] . we assume that the hidden charm mesons Z(4050) and .2^(4250) are vector 
tetraquark states, and study their masses using the QCD sum rules. The numerical 
results indicate that the masses of the vector hidden charm tetraquark states are about 
Mz = (5.12 ± 0.15) GeV or Mz = (5.16 ± 0.16) GeV, which are inconsistent with the 
experimental data and also much larger than the predictions of the constituent diquark 
model [I2l[l3l[l4l[l5]. 

The diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vector 
6*7^75, axial vector C7^ and tensor Ca^^. The structures 6*7^ and Ca^^ are symmetric, 
the structures C75, C and 6*7^75 are antisymmetric. The attractive interactions of one- 
gluon exchange favor formation of the diquarks in color antitriplet 3c, flavor antitriplet 3j 
and spin singlet 1^ |17t ll8j. The scalar hidden charm tetraquark states may have smaller 
masses than the corresponding vector states. 

The mass is a fundamental parameter in describing a hadron, in order to identify 
the Z(4050) and Z(4250) as tetraquark states, we must prove that the masses of the 
corresponding tetraquark states lie in the region (4.1 — 4.3) GeV. Furthermore, whether or 
not there exist such hidden tetraquark configurations is of great importance itself, because 
it provides a new opportunity for a deeper understanding of the low energy QCD. 

In this article, we assume that there exists a scalar hidden charm tetraquark state in 
the vr"^Xci invariant mass distribution, and construct the C75 — C75 type current Ji(x) 
and C — C type current J2ix) (and their superposition J{x)) to interpolate it, 

Ji(x) = e^^^^™"nJ(x)C75Cfc(x)c„(x)75CdI(x), (1) 

J2(x) = e^^'=e^™"nJ(x)Ccfc(a;)c„(x)C(il(x), (2) 

J{x) = cosO Ji{x) + sinO J2{x) , (3) 

where the i, j, ■ ■ ■ , n are color indexes; then study its mass using the QCD sum rules 
[13 [20]. The hidden charm mesons X(3872), y(4260), ^(4350), y(4660), Z(4430) have 
also been studied with the QCD sum rules as the tetraquark or molecular states [211 [221 
[231 [24]. 

In the QCD sum rules, the operator product expansion is used to expand the time- 
ordered currents into a series of quark and gluon condensates which parameterize the 
long distance properties of the QCD vacuum. Based on the quark-hadron duality, we can 
obtain copious information about the hadronic parameters at the phenomenological side 
[I3[20]. 

The article is arranged as follows: we derive the QCD sum rules for the mass of the 



Z in section 2; in section 3, numerical results and discussions; section 4 is reserved for 
conclusion. 

2 QCD sum rules for the tetraquark state Z 

In the following, we write down the two-point correlation function n(p) in the QCD sum 
rules, 

n(p) = i /"d^xe^P-^(0|r{j(x)jt(0)}|0), (4) 

we choose the scalar current J(x) to interpolate the tetraquark state Z . 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operator J{x) into the correlation function n(p) to obtain the 
hadronic representation |191 [20] . After isolating the ground state contribution from the 
pole term of the Z, we get the following result, 

"(^> = mF^ + -- (^' 

where the pole residue (or coupling) \z is defined by 

\z = (0|J(0)|Z(p)). (6) 

In the following, we briefly outline the operator product expansion for the correlation 
function Il{p) in perturbative QCD. The calculations are performed at the large space-like 
momentum region p^ <^ 0. We write down the "full" propagators Sij{x) and Cij{x) of a 
massive quark in the presence of the vacuum condensates firstly [20] , 
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where (qgsO'Gq) = {qgsO'af3G°'^q) and i ""^^ ) = { "" °^ ), then contract the quark fields 

in the correlation function Il{p) with Wick theorem, and obtain the result: 

U{p) = ifijk^imn^i'j'k'^i'm'n' j tfxe^^-^ {cOS^OTr [-f^Ckk'{x)-f^C S]^,{x)C] X 

Tr[^^Cr^,^{-x)^^CSl^{-x)C]+sin'eTr[Ckkix)CSl,{x)C] x 

Tr [Cm'm{-x)CSl,^{-x)C] ] . (9) 



Substitute the full u, d and c quark propagators into the correlation function Il{p) 
and complete the integral in the coordinate space, then integrate over the variables in 
the momentum space, we can obtain the correlation function Il{p) at the level of the 
quark-gluon degrees of freedom. 

We carry out the operator product expansion to the vacuum condensates adding up to 
dimension-10 and take the assumption of vacuum saturation for the high dimension vac- 
uum condensates, they are always factorized to lower condensates with vacuum saturation 
in the QCD sum rules, factorization works well in large Nc limit. In calculation, we observe 
that the contributions from the gluon condensate are suppressed by large denominators 
and would not play any significant roles [271 ESI [29l [30l |3l] . Furthermore, we neglect the 
terms proportional to the rriu and m^ as their contributions are of minor importance. 

Once analytical results are obtained, then we can take the quark-hadron duality and 
perform Borel transform with respect to the variable P^ = —p^, finally we obtain the 
following sum rule: 
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t = cos29e [-1,1]. 

Differentiating the Eq.(lO) with respect to -^, then eliminate the pole residue Xz, we 



can obtain a sum rule for the mass of the Z, 

M| = i5i^!^EW!)^^. (12) 

JZ^Jspis)e ^ 

3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 it 0.01 GeV)'^, 
{qgsCrGq) = ml{qq), ml = (0.8 ± 0.2) GeV^ (^^) = (0.33 GeV)^, and nic = (1.35 ± 
0.10) GeV at the energy scale about ^u = 1 GeV [El llQl ES]. 

The Belle collaboration observed the resonance-like structures Z(4050) and Z(4250) in 
the TT'^Xci invariant mass distribution near 4.1 GeV in the exclusive B^ — > K~tt~^Xc1 decays 
[ij. If they are scalar tetraquark states, the central value of the threshold parameter can 
be taken as sq = (4.248 + 0.5)^ GeV « 23 GeV , where we tentatively choose the energy 
gap between the ground states and the first radial excited states to be 0.5 GeV. 

The present experimental knowledge about the phenomenological hadronic spectral 
densities of the tetraquark states is rather vague, whether or not there exist tetraquark 
states is not confirmed with confidence, and no knowledge about the high resonances; we 
can borrow some ideas from the baryon spectra |26j . 

For the octet baryons with the quantum numbers I{J^) = \{\ )) the mass of the 
proton (the ground state) is Mp = 938 MeV, and the mass of the first radial excited state 
iV(1440) (the Roper resonance) is M1440 = (1420 - 1470) MeV ps 1440 MeV [26]. For the 
decuplet baryons with the quantum numbers I{J ) = |(| ) , the mass of the A(1232) 
(the ground state) is M1232 = (1231 - 1233) MeV Ri 1232 MeV, and the mass of the first 
radial excited state A(1600) is Mieoo = (1550 - 1700) MeV « 1600 MeV [26]. The energy 
gap between the ground states and the first radial excited states can be tentatively taken 
as 0.5 GeV for the light fiavor baryons. 

In Ref.[14j. Maiani et al assume the Z(4430) is the 25 1"' hidden charm tetraquark 
state to take into account the decay mode Z(4430) — > ip{2S) + vr^, as the IS" 1"' hidden 
charm tetraquark states lie at the interval (3750 — 3880) MeV and have the decay mode 
X^{1S) — > V'(l'S') + TT"*" or rjc{lS) + p'^ in the constituent diquark model [12]. The energy 
gap between the ground state and the first radial excited state is estimated to be M^(25) ~ 
^ip{iS) ~ 0.59 GeV for the heavy tetraquark states. 

We take it for granted that the energy gap between the ground states and the first 
radial excited states is about 0.5 GeV, and use this value as a guide to determine the 
threshold parameter sq with the QCD sum rules. 

We explore whether or not there exist scalar tetraquark states which consist of a 
scalar (pseudoscalar) diquark-antidiquark pair at the energy interval (4.1 — 4.3) GeV, and 
choose the larger value sq = (4.248 + 0.5)^ GeV ~ 23 GeV rather than the smaller value 
So = (4.051 + 0.5)^ GeV^ ss 22 GeV^ to take into account all possible contributions from 
the ground states. 

In the conventional QCD sum rules |19i I20j . there are two criteria (pole dominance 
and convergence of the operator product expansion) for choosing the Borel parameter M ^ 
and threshold parameter sq q 



^For the tetraquark states consist of light flavors, if the perturbative terms have the main contribution 
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Figure 1: The contributions from different terms with variation of the parameter t. The 
A and B correspond to the threshold parameters sq = 23 GeV^ and 25 GeV^ respectively. 
The notations a, /?, 7, A and r correspond to the perturbative term, {qq) + {qgsCrGq) term, 
(^) + (^) [{qq) + {qgsCjGq) + {qq?] term, {qq? + {qq){q9s<jGq)+ {qgscrGq? term 
and perturbative +{qq) + {qggCrGq) term, respectively. Here we take M^ = 3GeV^ and 
the central values of other input parameters. 
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Figure 2: The pole residue (or coupling) with variation of the parameter t. The notation 
a, /3 and 7 correspond to the threshold parameters sq = 23GeV^, 24GeV^ and 25GeV^, 
respectively. Here we take M"^ = 3 GeV and the central values of other input parameters.. 



The contributions from different terms with variation of the parameter t in the op- 
erator product expansion are shown in Fig.l. From the figure, we can see that the 
contributions from the term {qq) + {qggaGq) are negative at the interval — 1 < t < 0, 
which cancel out the contribution from the perturbative term greatly. The net contri- 
butions from the perturbative term-l-(^q') -|- {qggaGq) increase with variation of the t, 
and reach the largest value at t = 1. The contributions from the gluon condensates 
( ""^ ) + (^ — ) [{qq) + {QgsO'Gq) + {qq)'^] are very small and decrease with the param- 
eter t monotonously; the contributions from the high dimension condensates (qq)'^ + 
{qq){qgs(^Gq)+ {qgguGq)'^ also decrease with the parameter t monotonously. In other 
words, the operator product expansion converges more quickly for larger t at the interval 
t G [—1, 1], we can choose the value t = 1. 

On the other hand, the coupling of the interpolating current J(x) to the tetraquark 
state becomes stronger with larger t, see Fig. 2. It is reasonable to take the interpolating 
current with the strongest coupling to the tetraquark state. 

In Figs. 3-4, we plot the contributions from different terms in the operator product 
expansion. The contribution from the term ( "" ) is tiny and can be safely neglected. 
The contributions from the terms involving the gluon condensates are less than 8% even at 
very small Borel parameter M^, the gluon condensate plays a minor important role. The 



(in the conventional QCD sum rules, the perturbative terms always have the main contribution), we can 
approximate the spectral density with the perturbative term (where the A are some numerical coefficients) 

roo roo 

BuUr^A s*e~^ds = AM^° t^e'^dt, (13) 

Jo Jo 

then take the pole dominance condition, 



lo^t^e-^dt 



> 50% , (14) 



and obtain the approximated relation. 



to = ^>4.r. (15) 



The superpositions of different interpolating currents can only change the contributions from differ- 
ent terms in the operator product expansion, and improve convergence, they cannot change the leading 
behavior of the spectral density p{s) oc s* of the perturbative term [32j . 

This relation is difficult to satisfy for the light flavor tetraquark states [271 1281 [29l 1301 131) , if we take the 
Borel parameter has the typical value M^ = IGeV^, so > 4.7 GeV'^, the threshold parameter is too large 
for the light tetraquark state candidates /o(980), ao(980), etc. 

The hidden charm (or bottom) tetraquark states and open bottom tetraquark states may satisfy the 
relation, as they always have larger Borel parameter M^ and threshold parameter so [211 1221 1231 I24[ . 
Their spectral densities have the form p{s) = Cis* + C2S"^ + Css^ + • • • , where the d are coefficients, and 
exhibit the same leading behavior p(s) oc s^ as the light flavor tetraquark sates. If we take M^ = 1 GeV^, 
So > 4.7 GeV^, the threshold parameter so is too low for the hidden charm or open bottom tetraquark 
states, there is a large room for choosing larger threshold parameter to take into account the ground state 
contribution. We draw the conclusion that the hidden charm (bottom) tetraquark states and open bottom 
tetraquark states have possibility to satisfy the pole dominance condition. 

In this article, the vacuum condensate of the highest dimension {qggaGq)^ serve as a criterion for choosing 
the Borel parameter M^. At the value M^i„ > 2.2 GeV^, its contribution is less than 10% (see Fig. 3), we 
expect the operator product expansion is convergent. The relation in Eq.(15) indicates so > 10.5 GeV^, 
if we take a large Borel parameter M'^ > 2M^i„, then sq > 21GeV^, our phenomenological estimation 
So ~ 23 GeV^ is reasonable. 



vacuum condensate of the highest dimension {qgsCrGq)'^ serve as a criterion for choosing 
the Borel parameter M^. At the value M^-^ > 2.2 GeV , its contribution is less than 
10%, we expect the operator product expansion is convergent. 

The contributions from the vacuum condensates of high dimension {qq)"^ + {qq) {qgsCyGq) 
vary with the threshold parameter sq remarkably and serve as a criterion for choosing the 
threshold parameter sq. At the value sq > 23 GeV , their contributions are less than (or 
equal) 10% (see Fig.3-A), we expect the operator product expansion is convergent. The 
contributions from the vacuum condensates {qq)"^ + {qq){qgs(^Gq) + {qgsCrGq)'^ are less than 
13.5% at the values M^ > 2.2 GeV and sq ^ 23 GeV . The contributions from the vacuum 
condensates {qqf + {qq){qgs(jGq)+ {qgsaGqf+ (^) + (^) [{qq) + (^crGg) + {qqf] 
are less than 18%, the main contributions come from the perturbative term +{qq) +{qgs(jGq), 
see Fig. 4. The operator product expansion is convergent at the values Af^^^ > 2.2 GeV^ 
and So > 23GeV2. 

In Fig. 5, we plot the contribution from the pole term with variation of the threshold 
parameter sq. For the central values of the input parameters (except for t = 1), the 
contribution from the pole term is larger than 50% at the values M^„^. < 3.2 GeV and 
so> 23GeV2. 

In this article, the threshold parameter and the Borel parameter are taken as sq = (24ib 
1) GeV and M^ = (2.2 — 3.2) GeV respectively, the contribution from the pole term is 
about (51 — 88)%, the two criteria of the QCD sum rules are full filled JHHO]. We can take 
smaller Borel parameter and threshold parameter to satisfy the two criteria of the QCD 
sum rules marginally, however, the Borel window is rather small, M^^^ ~^min < -*- GeV^. 

Taking into account all uncertainties of the input parameters, finally we obtain the 
values of the mass and pole reside of the Z, which are shown in Figs. 4-5. From the 
figures, we can see that at the value M^ < 2.6 GeV^, the mass and the pole residue change 
remarkably with variation of the Borel parameter, we take the value M^ = (2.6—3.2) GeV , 
and obtain 

Mz = (4.36 ±0.18) GeV, 
\z = (3.38 ± 0.65) X 10^2 GeV^ (16) 

The meson Z(4250) may be a scalar tetraquark state, other possibilities, such as a hadro- 
charmonium resonances and a D^D^ + D^D^ molecular states are not excluded. 

The Z(4250) lie about (0.5 — 0.6) GeV above the vr+Xci threshold, if it is a tetraquark 
state, the decay Z -^ t^~^Xc1 can take place with the OZI super-allowed "fall-apart" mech- 
anism, which can take into account the large total width naturally; on the other hand, 
if it is a D^D^ + D^D\ molecular state, the decay can occur through the final-state re- 
scattering effects, Z — > D^D^ + D^D^ -^ t^'^XcI: and the corresponding width may be 
narrow, we have to search for other decay channels to accommodate the large total width. 

The typical decay mode Z -^ D^D^ is kinematically allowed, we can determine the 
spins of the Z(4250) with the angular distribution of the final state D^D^. If the decay 
Z -^ D^D^ is not observed (or the width is rather narrow), the Z(4250) may be a hadro- 
charmonium resonance (bound state of a relatively compact charmonium (xci) inside a 
light hadron (Tr"*") having a larger spatial size) [33]. The decay Z — > ir'^Xci occurs with the 
"fall-apart" mechanism and the width is large; while the decay Z -^ D^D^ takes place 
through the final-state re-scattering effects {Z — > tt^XcI -^ D^D^) and the width may be 
narrow. 
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Figure 3: The contributions from the vacuum condensates with variation of the Borel 
parameter M^. The ^4, i?, C, D, E and F correspond to the contributions from the 
(gg)2 +{qq){qgs(TGq) term, (qg^aGq)'^ term, (qq)'^ +{qq){qgs(^Gq) + (qgsCrGq)'^ term, 
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a, /3, 7, A, r and S, correspond to the threshold parameters sq = 20GeV^, 21GeV^, 
22 GeV^, 23 GeV^, 24 GeV^ and 25 GeV^, respectively. Here we take t = 1 and the central 
values of other input parameters. 
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4 Conclusion 

In this article, we assume that there exists a scalar hidden charm tetraquark state in 
the TT'^Xci invariant mass distribution, and study its mass using the QCD sum rules. 
The numerical result indicates that the mass is about Mz = (4.36 it 0.18) GeV, which 
is consistent with the experimental data. The hidden charm meson Z(4250) may be 
a tetraquark state. Other possibilities, such as a hadro-charmonium resonance and a 
DfD^ + D^D^ molecular state are not excluded; more experimental data are still needed 
to identify it. 
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